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■ Abstract 

We study the chiral loop corrections to the generalized Gerasimov-Drell-Hearn sum rule of the nucleon for 
finite photon virtuality in the framework of a Lorentz-invariant formulation of baryon chiral perturbation 
theory. We perform a complete one-loop calculation and obtain significant differences to previously found 
results based on the heavy baryon approach for the proton and neutron spin-dependent forward Compton 
amplitudes. 
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1. The spin structure of the nucleon is under active theoretical and experimental investigation. Of particular 
interest is the transition from the pcrturbative regime at high momentum transfer to the non-perturbative 
low-energy region. Interestingly, systematic and controlled theoretical calculations are only available for very 
small or very large momentum transfer while the intermediate region is accessible to resonance models or can 
be investigated using dispersion relations. Of special interest is the generalized Gerasimov-Drell-Hearn (GDH) 
sum rule [1] for virtual photons. It was pointed out in [2] that one can systematically investigate its momentum 
evolution at low energies in the framework of chiral perturbation theory (CHPT). In a series of papers, Ji et al. 
[3, 4] proposed a formalism that is better suited for the direct connection with the spin-dependent structure 
functions measured in deep inelastic scattering. A complete one-loop fourth order calculation was presented for 
the forward spin-dependent Compton amplitudes and it was shown that the Q 2 -variation was much stronger 
than expected from naive dimensional analysis. While the proton and the neutron spin-dependent amplitudes 
receive important resonance contributions that can not be calculated explicitly in chiral perturbation theory # 4 
it was pointed out by Burkert [6] that in the proton-neutron difference these resonance contaminations largely 
drop out and also that the Q 2 -variation is much more modest. The same happens at large Q 2 for the Bjorken 
sum rule in deep inelastic scattering which is free of certain assumptions entering the Ellis-Jaffe sum rules for 
protons and neutrons. Based on the results of [4], it was argued in [6] that the chiral prediction should only 
be taken seriously for Q 2 < 0.1 GeV 2 and that a continuous matching between the CHPT prediction and the 
pQCD evolution of the Bjorken sum rule to order a 3 can be questioned. Such type of smooth transition from 
the perturbative to the non-perturbative regime had been suggested from the study of the sum rules for nucleon 
polarizabilities at finite virtuality [7]. Note also that a vigorous experimental effort is underway at Jefferson 
Laboratory to determine the spin structure of the proton and the neutron. In this letter, we report about a 
novel analysis of the generalized Gerasimov-Drell-Hearn sum rule for low photon virtualities. We employ the 
Lorcntz-invariant formulation of baryon chiral perturbation theory proposed in [8], which in contrast to the 
heavy fermion approach automatically fulfills all strictures from analyticity (see also the earlier work in [9]). It 
furthermore has the advantage of resumming recoil corrections to all orders which in certain cases can improve 
the convergence dramatically, see e.g. the discussion of the neutron electric form factor in [10]. 

2. Next, we briefly summarize the formalism necessary for the discussion of the photon- nucleon Compton 
amplitude and the related sum rules. It is common to express the spin amplitude of forward doubly virtual 
Compton scattering in terms of two structure functions, called S\{v, Q 2 ) and S^^, Q 2 ), via 

T^\p, q, s) = -i e^q a {ro^Sifa Q 2 ) + [p ■ q Sp - S ■ q pp] } , (1) 

where s M denotes the spin-polarization four-vector of the nucleon, m is the nucleon mass, e^ ual3 the totally 
antisymmetric Levi-Civita tensor with e 0123 = 1, Q 2 > (the negative of) the photon virtuality and v = p-q/m 
the energy transfer. Note that while S\{v, Q 2 ) is even under crossing v <-> —v, S^^, Q 2 ) is °dd. In what follows, 
we will be concerned with the amplitudes S\(0, Q 2 ) = 5^(0, Q 2 ) — Sf(0, Q 2 ), i.e. the Compton amplitudes with 
the contribution from the elastic intermediate state subtracted. The generalized sum rule proposed in [2] utilizes 
the imaginary part of the photon-nucleon spin-flip amplitude fe(v, Q 2 ), which can be expressed in terms of the 
combination S\(0, Q 2 ) — (Q 2 /m)dS2(v, Q 2 )jdv\ v=a # 5 . We concentrate here on the sum rule [4] 

/J V Gl( "'° 2) = ^s^Q 2 ) = 4? HQ 2 ) , (2) 

with G\(v,Q 2 ) a spin-dependent structure function, G\(v,Q 2 ) — g\(Q 2 /2mv, Q 2 )/(mv) and g\ the dimension- 
less scaling function. Furthermore, denotes the inelastic threshold. The subtracted Compton amplitude 
Si(v, Q 2 ) is normalized as Si (0,0) = —e 2 K 2 /m 2 , with k the anomalous magnetic moment of the spin-1/2 par- 
ticle and e 2 /47r = 1/137.036 the fine structure constant. This formula can be converted into the first moment 
of the sum rule via 

r 1 (Q 2 ) = ^7(Q 2 ). (3) 

# 4 An exception to this is the contribution from the A(1232) which can be systematically analyzed in the framework of the small 
scale expansion [5]. However, at present a relativistic formulation for the A based on a systematic power counting has not been 
worked out. 

# 5 Note the different convention as compared to Ref. [3]. 
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All formulas given so far hold for the proton and the neutron separately and the corresponding quantities will 
be indicated by a superscript p or n. At small photon virtualities, the subtracted Compton amplitudes §i(0, Q 2 ) 
(i = 1,2) can be calculated in chiral perturbation theory, as will be discussed next. For more details on the 
formalism and earlier calculations we refer to [3] . 

3. We briefly discuss the effective chiral Lagrangian underlying our calculation. Its generic form consists of a 
string of terms with increasing chiral dimension, 

£ eff =4 1 U^+4 3 U4 4 U4 2 i+4 4 »+- • (4) 

The superscript denotes the power in the genuine small parameter q (denoting pion masses and/or external 
momenta). A complete one-loop (fourth order) calculation must include all tree level graphs with insertions 

(2) 

from all terms given in Eq. (4) and loop graphs with at most one insertion from C\. ^ . The complete Lagrangian 
to this order is given in [11]. We note that for the case under consideration the only appearing dimension two 
LECs, called c 6 and cj [12], can be fixed from the anomalous magnetic moment of the proton and of the neutron, 
respectively. 

Baryon chiral perturbation theory is complicated by the fact that the nucleon mass does not vanish in the chiral 
limit and thus introduces a new mass scale apart from the one set by the quark masses. Therefore, any power 
of the quark masses can be generated by chiral loops in the nucleon (baryon) case, spoiling the one-to-one 
correspondence between the loop expansion and the one in the small parameter q. One method to overcome 
this is the heavy mass expansion (called HBCHPT) where the nucleon mass is transformed from the propagator 
into a string of vertices with increasing powers of 1/m. Then, a consistent power counting emerges. However, 
this method has the disadvantage that certain types of diagrams are at odds with strictures from analyticity. 
The best example is the so-called triangle graph, which enters e.g. the scalar form factor or the isovector 
electromagnetic form factors of the nucleon. In a fully relativistic treatment, such constraints from analyticity 
are automatically fulfilled. It was recently argued in [9] that relativistic one-loop integrals can be separated into 
"soft"' and "hard" parts. While for the former the power counting as in HBCHPT applies, the contributions 
from the latter can be absorbed in certain low-energy constants (LECs). In this way, one can combine the 
advantages of both methods. A more formal and rigorous implementation of such a program was given in [8]. 
The underlying method is called "infrared regularization" . Any dimcnsionally regularized one-loop integral H is 
split into an infrared singular and a regular part by a particular choice of Feynman parameterization. Consider 
first the regular part, called R. If one chirally expands these terms, one generates polynomials in momenta and 
quark masses. Consequently, to any order, R can be absorbed in the LECs of the effective Lagrangian. On the 
other hand, the infrared (IR) singular part I has the same analytical properties as the full integral H in the 
low-energy region and its chiral expansion leads to the non-trivial momentum and quark-mass dependences of 
CHPT, like e.g. the chiral logs or fractional powers of the quark masses. For a typical one-loop integral (like 
e.g. the nucleon self-energy) this splitting can be achieved in the following way 

f d d k 1 f 1 f d d k 1 f / >oc _ [°°\ f d d k 1 

H J (2ny AB~J Z J (2ny [(1 - z)A + zB} 2 \J J 1 J Z J (2tt)<* [(1 - z)A + zB} 2 + K ' 

with A = M 2 — k 2 — ie, B = m 2 — (p — k) 2 —ie, e — > + andM„ the pion mass. Any general one-loop diagram with 
arbitrary many insertions from external sources can be brought into this form by combining the propagators 
to a single pion and a single nucleon propagator. It was also shown that this procedure leads to a unique, 
i.e. process-independent result, in accordance with the chiral Ward identities of QCD [8]. Consequently, the 
transition from any one-loop graph H to its IR singular piece / defines a symmetry-preserving regularization. 

Within this approach, we have calculated S[ P i n) (0,Q 2 ). The pertinent one-loop diagrams are given in [4], 
however, in the approach used here one does not need to consider dimension two insertions from the kinetic 
energy on the nucleon lines since this is done automatically in a Lorentz-invariant formulation. In fact, all such 
diagrams are counting as third order in the approach used here, the only genuine fourth order graphs are the 
ones with one insertion ~ cqj (anomalous magnetic moment). Therefore, S\(u, Q 2 ) is already non-vanishing at 
third order, quite in contrast to the heavy baryon calculation. As a check, we remark that in the limit of a very 
large nucleon mass, one recovers the earlier HBCHPT results of [2] and [4]. However, as already found in the 
relativistic calculation of the slope of the generalized GDH sum rule in [2] , one can not give the results in closed 
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Figure 1: Chiral loop contribution to the structure function 5i(0, Q 2 ) with the elastic contribution subtracted. 
The solid (dot-dashed) line gives the result of the present calculation to order q 4 (q 3 ) in comparison to the 
heavy baryon result of [4] (dotted line). Left (right) panel: Proton (neutron). 



analytical form. Wc refrain from giving the expressions based on scalar loop functions here but rather refer 
to the upcoming paper [13]. Of course, there are also non-negligible resonance contributions to the Si(v,Q 2 ). 
More precisely, the effects from the A largely cancel in Si(v, Q 2 ) [4] but are expected to be more important i 



m 



Si{y, Q 2 )- There are additional (smaller) vector meson and higher baryon resonance contributions. Therefore 
we mostly concentrate here on S\(y, Q 2 ) and also on the neutron-proton difference in which most resonance 
effects cancel, see e.g. [14, 15]. 

4. We now present our results. First, we must specify the parameters. We use §a = 1.267, F„ = 93MeV, 
M v = 139.57 MeV, to = 0.9389 MeV, k v = 3.706 and k s = -0.120. The latter two quantities include the LECs 
c 6 and c 7 . We point out again that there are no unknown parameters and thus we obtain genuine one-loop 
predictions. In Fig. 1 we show the results for 5f' n (0, Q 2 ) in comparison to the heavy baryon result of [4]. Some 
remarks are in order. First, we note that in the IR approach, one has third and fourth order contributions 
with the latter becoming sizeable for Q 2 > 0.2 (0.3)GeV 2 for the proton (neutron) #6 . Second, the modulus 
of the structure function is bigger as compared to the HBCHPT result. Similar remarks apply to the results 
for 5*2 '"(0, Q 2 ) not shown here. These sizeable differences between our results and the HBCHPT ones can be 
traced back to the fact that the 1/m expansion for these structure functions is very slowly converging, and thus 
the method used here appears to be more appropriate. From these results one deduces the following slopes for 
the proton and neutron structure function S\{Q 2 ) (from here on, we do not display the argument v = any 
longer) : 

d ^ (g2) = (-0.4 + 1.2) GeV- 4 -0.8 GeV- 4 , 

Q 2 =0 



dQ 2 
dS?{Q 2 ) 



dQ 2 



= (-1.7 + 1.9) GeV~ 4 = 0.2 GeV" 4 , (7) 

Q 2 =0 



where the first (second) term in the round brackets refers to the third (fourth) order contribution. These results 
are compatible with expectations from naive dimensional analysis. For comparison, the corresponding HBCHPT 
results are dSf/dQ 2 = 2.9GcV~ 4 and dS'?/dQ 2 = 2.4GeV~ 4 , which are quite different. The large slope in the 
heavy baryon case is due to a variety of effects. First, the large value of the isovector anomalous magnetic 
moment, k v = 3.71, appearing in dSi/dQ 2 oc [1 + 3k v + 2(1 + 3k s )t 3 ] leads to a prefactor for the proton of 14 
instead of the expected order of one based on naive dimensional analysis. Second, as as been spelled out clearly 
in [16], in the spin sector the natural expansion parameter seems to be -KM^/m rather than M^/m (for the odd 
powers), which can lead to sizeable corrections. Consequently, the 1/m expansion converges slowly, and thus 

# 6 We remark that for certain tree graphs at v = the heavy baryon propagator formally scales as 0(p~ 2 ) so that one could 
consider the contribution generated from such diagrams also as third order in HBCHPT, for details see [13] 
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Figure 2: Chiral loop contribution to the first moment Yi(Q 2 ). The solid (dot-dashed) line gives the result 
of the present calculation to order q 3 (q 4 ) in comparison to the heavy baryon result of [4] (dotted line). Left 
(right) panel: Proton (neutron). 



the numbers for these slopes found here should be considered more reliable. Note also that for these slopes 
the relative size of the fourth order corrections comes out very different for the proton and the neutron, but 
it should be stressed that the meaningful numbers are the full one-loop (third and fourth order) results. For 
orientation, we quote the result of the dispersive analysis of [17], which besides pion loop effects also includes 
resonance contributions, it is dS^/dQ 2 = 1.8 GeV -4 . 

In Fig. 2 we show the first moments T\{Q 2 ) and T"(Q 2 ) in comparison to the heavy baryon results. For the 
range of photon virtualitics below Q 2 = 0.2 (0.3) GeV 2 , the fourth order corrections are not dramatically large 
for the proton (neutron). Note that while T^(Q 2 ) crosses zero at Q 2 ~ 0.33 GeV 2 , T™(Q 2 ) stays negative for all 
values of Q 2 in the range of virtualities considered here. We stress again that for a direct comparison with the 
data, some resonance contributions have to be added, this will be done and discussed in the upcoming work 
[13]. However, one does not expect these resonance contributions to be overly large. Note that we refrain from 
comparing these curves with the preliminary data from Jefferson Laboratory because these data are not yet 
final and can therefore not be shown. This comparison will be done in [13]. 

In Fig. 3 we show the result for T\~ n {Q 2 ) in comparison to the HBCHPT result [6]. This should be considered 
the main result of this paper. For Q 2 < 0.1 GeV 2 , the fourth order IR and the heavy baryon result are very 
similar. For larger momentum transfer, however, Y\~ n (Q 2 ) grows faster with increasing momentum compared 
to the heavy baryon calculation. This is due to the fact that the chiral corrections to the individual T\(Q 2 ) and 
r™(Q 2 ) are quite different (similar) in the IR (HB) approach (for Q 2 > 0.1 GeV 2 ). The corresponding slope is 
of course given by the GDH sum rule, 



-2 ,„2 



drr i (Q 2 ) _«* «i = 063 GeV - 2j (8) 



dQ 2 q 2 =o 8m 2 
as can also been seen in the figure since all curves merge to one around Q 2 ~ 0. 

It is also interesting to consider the slopes corresponding to the combination SuiQ 2 ) = <Si(0, Q 2 ) — (Q 2 /m) 
x(dS i 2(i / , Q 2 ) / dv)\ u= Q for the proton and the neutron case separately, we find 

dSUQ 2 ) 



dQ 2 
dSUQ 2 ) 



dQ 2 



= (3.2 + 1.8) GeV = 5.0 GeV , (9) 

Q 2 =0 

= (-0.3 - 1.8) GeV" 4 = -2.1 GcV~ 4 , (10) 

Q 2 =0 

which should be compared with the heavy baryon results [2, 4], dSf 2 (Q 2 )/dQ 2 = 3.7(1 — 2.6) GeV~ 4 = 
-6.0 GeV" 4 and dS[ l 2 (Q 2 )/dQ 2 = 3.7(1 - 2.1) GeV" 4 = -4.3 GeV" 4 . We remark that the seemingly over- 
whelming fourth order correction for the neutron is simply a reflection of the accidentally suppressed third order 
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result. The 50% correction found for the proton is not uncommon for the nucleon spin sector. The larger cor- 
rections in HBCHPT are simply a reflection of the various effects discussed already before, In the IR approach 
used here, higher order 1/m corrections (in the HBCHPT counting) tend to diminish such large corrections. 
We point out that such effects also appear in other observables, like e.g. in pion photoproduction off nucleons, 
but in that case their effect is much less dramatic because they are accompanied by small or normal prefactors 
(less than or of order one) or do not play a significant role. Such large corrections are only appearing in the 
spin sector and their origin is essentially understood. Note further that the slope of the integrals I P: n{0) defined 
in [2] differ from the conventions used here by an overall factor of n/2. For this combination of the structure 
functions, one has to include the A (and other resonance) contribution for a direct comparison with the data, 
see e.g. [2]. 

5. In this paper, we have considered the one-loop representation of the structure functions of doubly virtual 
Compton scattering at low photon virtualities, Si(0, Q 2 ), i = 1,2 (with the elastic intermediate state contribu- 
tion subtracted). We have performed the analysis in the framework of infrared regularized relativistic baryon 
chiral perturbation theory and compared the results with the existing ones obtained in the heavy mass formal- 
ism. We have found that there are significant differences for the individual structure functions Si(0,Q 2 ) and 
5f(0, Q 2 ), i.e. these exhibit a less dramatic (stronger) Q 2 -evolution for the proton (neutron) than in HBCHPT 
and show reasonable convergence for small photon virtualities. We have also pointed out that the 1/m expan- 
sion underlying the heavy baryon approach is only slowly converging for these structure functions and thus 
the resummation of recoil effects inherent to the IR formalism is to be considered important. Clearly, these 
predictions should only be compared to data when the missing resonance contributions are included. However, 
in the first moment difference r^™(<5 2 ), contributions from spin-3/2 resonances like the delta drop out (and 
most other resonance contributions, too), so that the predictions shown in Fig. 3 should be considered the main 
result of this paper. As it turns out, the growth of T\~ n (Q 2 ) with increasing Q 2 is stronger to what was found 
previously in the heavy mass formulation. It would be interesting to compare this prediction with the data soon 
becoming available from Jefferson Laboratory^ 7 . It is important to systematically add resonance contributions 
to these chiral predictions for the proton and the neutron separately to allow for a direct comparison with the 
data. Work along these lines is underway [13]. 

# 7 First preliminary data from the CLAS collaboration which start at Q 2 = 0.2 GeV 2 seem to be consistent with our and the 
heavy baryon result [18]. 
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Note added 

After submission of this manuscript, new experimental information on the integral I n (Q 2 ) for Q 2 = 0.1 — 
0.9 GeV 2 appeared (M. Amarian ct al. [The Jefferson Lab E94010 Coll.], arXiv:nucl-ex/0205020). The point 
at lowest photon virtuality is in good agreement with the prediction given here, but the prediction at Q 2 = 
0.26 GeV 2 is much more negative than the data. 
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